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Abstract: We derive a formula for the entanglement entropy of squeezed states on a lattice 
in terms of the complex structure J. The analysis involves the identification of squeezed states 
with group-theoretical coherent states of the symplectic group and the relation between the 
coset Sp{2N,Mf)/Isot{Jo) and the space of complex structures. We present two applications 
of the new formula: (i) we derive the area law for the ground state of a scalar field on a 
generic lattice in the limit of small speed of sound, (ii) we compute the rate of growth of 
the entanglement entropy in the presence of an instability and show that it is asymptotically 
bounded from above by the Kolmogorov-Sinai rate. 
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1 Introduction 

Squeezed states play a prominent role in the description of a wide range of phenomena in¬ 
cluding non-classical states of light in quantum optics [1], the emission of radiation from 
evaporating black holes [2] and the amplification of quantum fluctuations into macroscopic 
perturbations during cosmic inflation [3-7]. In this paper we consider squeezed states for a 
system of N bosonic degrees of freedom. The system can represent for instance the lattice 
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discretization of a scalar field. Given a set of creation and annihilation operators a*, aj lo¬ 
calized at each site of the lattice, we have a notion of Fock vacuum and a family of squeezed 
vacua defined as the set of states obtained by applying the squeeze operator exp(^Y^alaj) to 
the Fock vacuum [1]. In fact each squeezed vacuum can be understood as the Fock vacuum 
associated to some set a*, oj of creation and annihilation operators obtained from the first 
set via a Bogolyubov transformation [8]. The mathematical structure behind this ambiguity 
is the choice of a complex structure J: a squeezed vacuum | J) is labeled by an element of 
the space of positive compatible complex structures of the system. We derive a new formula 
for the entanglement entropy of a squeezed state: given a decomposition of the Hilbert space 
of the system in a tensor product 'Hy = "Ha ® the entanglement entropy 5^4(1 J)) of a 
squeezed vacuum restricted to the subsystem T-La is 


SA{\J))=tT 


\ - PaUPa , 
--log 


Pa'iJPa 
2 


where Pa is the projector to the phase space A of the subsystem and the trace is over N 
dimensional matrices. This formula has a clear interpretation: there is entanglement between 
A and B only if the projection Pa JPa of the complex structure is not a complex structure 
for the subsystem. In the derivation we use in a crucial way the map between the space 
of complex structures and the coset Sp{2N,'R)/Isot{Jo) of symplectic transformation that 
squeeze Jq. This map allows us to formulate squeezed vacua as group-theoretical coherent 
states for the symplectic group [9] and to determine their entanglement Hamiltonian. 


The ground state of a free scalar field discretized on a lattice provides a simple example 
of squeezed vacuum. Its entanglement entropy has been computed in the classical works 
[10-12] and shown to result in an area law [13]. More recently upper and lower bounds on 
the area scaling have been determined using results on block banded matrices [14-17]. As an 
application of our formula, in Sec. 5.2 we determine the entanglement entropy of the ground 
state of a free scalar field Hamiltonian discretized on a generic lattice and find the expression 
Sa{\J)) = T(logD Area(aA) +0(0^) in the limit of small speed of sound c. 

The states we are interested in form a larger class than just ground states of quadratic 
Hamiltonians with local interactions. Squeezed vacua can be equivalently characterized in 
terms of their 2-point correlation function: in the Schrodinger representation they are complex 
Gaussian functions of N variables. Gaussian states play a central role in quantum information 
theory as they constitute a versatile resource for quantum communication protocols [18]. 
Their entanglement entropy was derived in [19] and later formulated in terms of symplectic 
invariants [20]. The analysis of this paper largely relies on the methods developed in this 
setting (see [21-24] for recent reviews). The focus here however is not on the covariance matrix 
but on the complex structure J. Our motivation comes from the possible field-theoretical 
applications of our formula. 

The complex structure plays a central role in the definition of the notion of particles 
in quantum field theory in curved space [25-28] and in the description of unitarily inequiv¬ 
alent representations of the canonical commutation relations [8, 29]. In the case of black 
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hole evaporation the paradigmatic setting involves an initial state prepared in the in-vacuum 

1 Jin), its evolution in a black hole space-time and its analysis in terms of out-particles defined 
by the complex structure Jout, [2, 25-28]. At late times the evolution of the state |Jin) is 
a squeezed vacuum with respect to the complex structure Jout- Its entanglement entropy 
and thermal properties have mostly been studied using conformal field theory techniques in a 
two-dimensional approximation [30-32]. In the case of inflationary cosmology [3, 4], an initial 
state \ Jbd) evolves in a quasi de Sitter spacetime. The roll-over transition results in an insta¬ 
bility in the evolution Hamiltonian and at a cosmic time r] before the end of inflation the state 
evolves in a highly squeezed vacuum | J^), [3? -6]. In this paper we do not address directly 
the question of the definition and computation of the entanglement entropy of such squeezed 
states in quantum field theory: we restrict attention to systems with finitely many degrees 
of freedom, effectively introducing a lattice cut-off of the continuum theory. In this setting 
we study the evolution of a squeezed state in the presence of an instability. The instability 
is modeled using a quadratic lattice Hamiltonian that is unbounded from below. In the long 
time limit we find the remarkably simple result reported in Sec. 5.3: the rate of growth of the 
entanglement entropy is asymptotically bounded from above by the Kolmogorov-Sinai rate 
of the associated classical system [33] and is only mildly dependent on the initial squeezed 
state and the choice of subsystem A. 

We begin by defining the structures involved in the derivation of our main formula for 
the entanglement entropy 5 a (|J)) of a squeezed vacuum |J). The main formula is derived 
in Sec. 4 and used in section Sec. 5 to study the area law and the time-evolution of the 
entanglement entropy in the presence of instabilities. 

2 Symplectic group and N bosonic degrees of freedom 

In this section we summarize some basic results about the symplectic group and its unitary 
representation on the Hilbert space of a system of N bosonic degrees of freedom [34-37]. 

2.1 Symplectic vector space (H, H) and the space of complex structures 

Consider a linear vector space V = of finite dimension 2N. We denote the elements 
of V and Wa the elements of its linear dual V''. A symplectic structure on H is a bilinear map 
^ab ■ y X H — )• M that is antisymmetric 

^ab V°'U^ = -Qab U°‘v’^ ( 2 . 1 ) 


and non-degenerate, i.e. 


= 0 ^ u“ = 0. (2.2) 

The inverse is denoted : Wx )- M and satisfies the relation where 

is the identity. The couple (H, ^ab) is a symplectic vector space of dimension 2N. 
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A linear subspace A C is a symplectic subspace of V if the intersection with its 
symplectic complement 


= {v'^ g i/| = 0, Vu'’ G A} (2.3) 

is null, i.e. A n A^ = {0}. Defining B = A^ we have the decomposition V = A Q) B. We 
denote Pa '■ V ^ V the projector from V to the subspace A. 

A linear symplectic transformation is a linear map : V ^ V that preserves the 
symplectic structure 

(2.4) 

or equivalently in terms of the inverse symplectic structure Linear sym¬ 

plectic transformations M^h form a group, the linear symplectic group Sp{2N,M.). 

On the symplectic vector space V we introduce a further structure, a linear complex 
structure, i.e. a linear map J“;, : V ^ V that squares to minus the identity, 

J\j\ = -5\. (2.5) 

The complex structure J“;, is said to be compatible to the symplectic structure VLab and 
positive if the following two conditions are satisfied: 

^abJ°'cJ^d = ^cd, ( 2 . 6 ) 

> 0. (2.7) 

The first equation tells us that the complex structure belongs to the linear symplectic group, 
G Sp{2N,'R), the second that together with the symplectic structure it defines a metric, 
i.e. a positive definite inner product gab : D x D —)■ M given by 

9ab = ^acJ'^b- ( 2 . 8 ) 

We are interested in the set of all complex structures on {V, D) that are compatible and 
positive. We denote this set 

J (D, D) = {I positive compatible complex structure on {V, D)} . (2.9) 

Note that the projection to A of a complex structure, i.e. PaJPa, in general does not induce 
a complex structure on A. 

2.2 Polar decomposition and the Squeeze{jQ) coset 

In the following we introduce a matrix notation for symplectic transformations and their 
polar decomposition. It is evident from its definition that this notation involves the choice of 
a reference complex structure Jq. 
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We denote a vector v^' € V simply by u = (f“). Similarly a linear map ; y —)• 1/ is 
denoted M = (M";,). The metric g^ab = ^ac Jifb can be understood as a linear map from the 
vector space to its dual, g^ab '■ V ^ V''. We use this map to lower indices, i.e. Ua = goabU^, 
and define the transpose = {goabu’^)-^ 

Note that the definition of the transpose involves the choice of a reference complex struc¬ 
ture Jo = (Jo^b)- For definiteness we choose a reference complex structure Jo associated to 
a decomposition of the symplectic vector space V = in symplectic subspaces Vi = 
so that V = (BiVi with i = 1,... ,N. The restriction flj of the symplectic structure n to 
each subspace Vi can be put in the canonical form with an appropriate choice of coordinates. 
Given a positive compatible complex structure jo in each Vi we define Jq as the direct sum 
over i. In particular we choose 

io = (J o), = = (i 0 )’ 

where 1 = (dj) is the N x N identity. It is immediate to check that the two conditions 
(2.6) and (2.7) are satisfied so that the complex structure is compatible with the symplectic 
structure and positive. With this choice of complex structure, the associated metric (2.8) 
is simply the Euclidean metric goab = ^acJo^b = ^ab ■ In particular in the matrix notation 
discussed above the transpose of a vector is simply = (<!„{, tt^). 

In matrix language, the complex structure Jq satisfies = —1 and Jo* = — Jo- The 
linear symplectic group is the group of 2N x 2N matrices satisfying MJq = Jq, i.e. 

Sp{2N, R) = {M G Mat(2iV, M) | MJq M* = Jo} . (2.11) 

From the definition it is immediate to prove that a symplectic matrix M is invertible and it 
has inverse^ 

M“* = -JoM*Jo. (2.12) 

The linear symplectic group Sp{2N,R) is non-compact. A maximally compact subgroup 
can be identified with the set of elements of Sp{2N, M) that leave the chosen complex structure 
Jo £ J{V,Vt) invariant under conjugation, i.e. the isotropy group Isot{Jo) defined as 

Isot{Jo) = {M G Sp{2N,R) I MJqM-^ = Jq] . (2.13) 

Using (2.12) we see that the elements of Isot^Jo) satisfy MM* = 1, i.e. they are orthogonal 
transformations in 0{2N,R) and 

Isot{jQ) = Sp{2N,R) n 0{2N,R). (2.14) 

^This notation is consistent with being the inverse of Blab- 

^It is less immediate to prove that detM = -|-1. An algebraic proof makes use of the notion of PfafBan. 
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As 0{2N,R) is a compact group, it follows that also the isotropy group is compact, in fact it 
is a maximally compact subgroup of Sp{2N,R). 

Symplectic transformations that ‘squeeze’ the complex structure Jq, i.e. that do not leave 
it invariant under conjugation, form the coset space 

Squeeze{Jo) = Sp{2N,'R)/Isot{jQ). (2-15) 

The polar decomposition allows us to write uniquely every symplectic matrix as a product of 
a matrix R G Isot{Jo) and a matrix T G Squeeze{Jo), 

M = TR, (2.16) 

where T = (M M^) 2 and R = It is clear that RR^ = 1 and therefore R G 

0{2N,R). The matrix T = {M M ^)2 is symmetric and positive, T G Sym_|_(2A', M) and 
provides a choice of representative in the squeeze coset 

Squeeze{Jo) ~ Sp{2N,R) n Sym+(2A^,M). (2.17) 

The Lie algebra sp(2A^, M) of the symplectic group and its decomposition in generators of 
Isot{Jo) and of Squeeze{Jo) is discussed in App. A. 

The space of positive compatible complex structures J' (17, fl) can naturally be identified 
with the Squeeze{Jo) coset, 

~ Squeeze{Jo). (2-18) 

In fact every J G can be written as J = MJqM~^ for some M G Sp{2N,R) and 

Squeeze{jQ) is the set of symplectic transformations that act non-trivially on Jq under con¬ 
jugation. 

2.3 Phase space and Poisson brackets 

The symplectic vector space (P, fl) can be understood as the phase space of a classical dy¬ 
namical system with N degrees of freedom. Consider coordinates on P = and smooth 
functions f,g G C'“(M^^). The Poisson bracket 

{f,9} = ^'^^dafdbg (2.19) 

endows the space of smooth functions with a Lie algebra structure. 

The space of constant or linear functions on P, i.e. functions of the form /(.^“) = 
c + forms a sub-algebra called the Heisenberg algebra f) 2 Af-i-i- Canonical variables are 
linear functions g* = p* = such that they satisfy the canonical Poisson brackets 

= 0, = 0 with i,j = In more compact form, by 

viewing the coordinates as functions on phase space, we can write the canonical Poisson 
brackets as 
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( 2 . 20 ) 


Linear canonical transformations are transformations —)• M°‘i, that preserve the canonical 
Poisson brackets, i.e. Therefore the matrix M“6 is an elements of the 

symplectic group 5p(2iV, M). 

The space of quadratic functions on V, i.e. functions of the form 

H = lKabee, ( 2 . 21 ) 

also forms a subalgebra: 

{Hi,H2} = H3 ( 2 . 22 ) 

where = Q^KiacK 2 bd This algebra is isomorphic to the Lie algebra sp(2A^, M) of 

the symplectic group (App. A). Quadratic functions on phase space are the generators of 
linear canonical transformations. 

2.4 Quantization, the Hilbert space Tiv, creation and annihilation operators 

Quantization of the symplectic vector space iV, O) results in a Hilbert space Hy of N bosonic 
degrees of freedom. The elementary observables in the quantum theory satisfy the canonical 
commutation relations 

= i^l“^ (2.23) 

The Hilbert space T-iy provides a regular representation of the Heisenberg algebra i)2N+i and 
can be realized in the Schrodinger representation as the space L^(M'^) of square-integrable 
functions depending on a maximally commuting set of linear observables, e.g. V’(9*)' For finite 
N, the Stone-von Neumann theorem states that irreducible representations of the canonical 
commutation relations on a Hilbert space are unitarily equivalent to the Schrodinger repre¬ 
sentation [38]. In the following we adopt the Fock representation defined in terms of creation 
and annihilation operators [8]. 

The reference complex structure Jq defined in Eq. (2.10) can be diagonalized in Vc = 
H (g) C = with a unitary transformation^ and put in the form 

N 

= (2.24) 

i=l 

where Ui and u* are orthonormal eigenvectors in with 

< = ^ - i ^i+N) ■ (2-25) 

The vectors span the eigenspace Vc"*" C Vc corresponding to the eigenvalue -|-i of Jq, and 
analogously u^* spans the eigenspace Vc~ with eigenvalue —i. Creation and annihilation 
operators associated to the complex structure Jq are defined as 

ai = uie, 4 = Uiae- (2.26) 

^Note that iJo is a Hermitian matrix. 
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Using the expression 1 = identity we can write the operator as 


N 

e = Y.(-‘^‘la, + uta\). ( 2 . 27 ) 

i=l 

Thanks to the compatibility of Jq to the symplectic structure, the creation and annihilation 
operators satisfy the canonical commutation relations 

[oj, a}-] = 6ij , [ui, aj] = 0 , [a\, aj] = 0. (2.28) 

The Hilbert space Tiy is defined algebraically by introducing the reference vacuum | Jq) defined 
as the state annihilated by all the Oj 

Oi|Jo) = 0 i = (2.29) 

and the orthonormal set of states 

N / "j" \ 

Im,... ,njv; Jo) = ^ ']\Jo) , (2.30) 

with Ui G Nq. The Hilbert space T-lv is the completion of the span of the Fock basis jn,; Jq) = 
|ni,..., UN', Jo)- A generic state |s) G % can be expanded on the basis (2.30) and is given by 

\s) = '^Cni\ni;Jo), (2.31) 

rii 

with c„. a complex square-summable sequence. 


The relation between the complex structure Jq and the state | Jq) can be made manifest 
by noticing that the matrix 

k — l-ko X—^ t / \ 

= (2.32) 

i=l 

is the projector from Vc = Vc"'' © Vc~ to Vc^. Therefore Eq. (2.29) can be written as 


k — i Jq 
2 


i\Jo) = o. 


(2.33) 


Note that a different choice of complex structure J va. J (U, H) corresponds to a different 
choice of state | J) that can be used as reference vacuum for the Fock basis. The state | J) is 
defined by 


1 




analogously to the vacuum | Jo) associated to the complex structure Jq. 


(2.34) 






2.5 Unitary representation of Sp{2N,'R.) on the Hilbert space Tiv 


The Hilbert space T-Ly carries a unitary representation of the symplectic group Sp(2A^, M),^ 
[36, 37] 

[/(M)^“U(M)"^ = (2.35) 

Using the canonical commutation relations it can be shown that a quadratic hermitian oper¬ 
ator with Kab £ Sym(2A^, M) generates the unitary transformation representing a 

symplectic transformation M with generator Q.°‘^Kcb, i-e. 

U{M) = exp (-Z with M“b = exp{Q^^K^b) ■ (2.36) 

Symplectic transformation of the creation and annihilation operators associated to the refer¬ 
ence complex structure Jq follow from Eq. (2.26) and are given by 

U{M) ai U{M)-^ = aj + a] (2.37) 

with the Bogolyubov coefficients given by = u\Muj and 'hij = ujMnt, ^ [8]. The set of 
quadratic hermitian operator that generate these unitary transformations can also 

be expressed in terms of creation and annihilation operators: we define the operators Eij, Fij 
and fJj 

Fij = “(ojoj + o.jO'i )) (2.38) 

Fij = ai aj , f}- = 4 • (2.39) 

Their commutation relations provide a representation of the symplectic algebra sp(2A^, M), 


1 

II 

(2.40) 

[Fij, Ffti] — 0 , [T)4 F^ij — 0 

(2.41) 

[F^ij^ Ej^i\ — Fii5j]^ -|- FjiSij^ 

(2.42) 

[Fl,E,i] = -Fl6ji-Fl6u 

(2.43) 


(2.44) 


'^The unitary operators U{M) provide in fact a projective representation of the symplectic group. A double 
covering of the group is required to make the representation single valued. This double covering is called the 
metaplectic group. An example of the two-valued nature of the representation is the following. Consider the 
element Mt = exp(tJo) of Sp{2N,'R). We have Mg = Mg-n = 1- On the other hand the unitary representation 
of Mt is given by U{Mt) — exp(it with the operator having a half-integer spectrum for N 

odd. As a result, while U{Mg) = 1 we can have U{M 2 tt)\s) = ±|s) on eigenstates of 

® This relation can be inverted so that the symplectic transformation M is written in terms of the Bo¬ 
golyubov coefficients as M = -g 
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From Eq. (2.40) we see that the Eij form a subalgebra: they are the generators of the unitary 
group U{N), i.e. 

exp(iQ;®'^ Eij) ai exp(— Eij) = Eij aj (2-45) 

with Eij a unitary matrix. The transformations generated by Eij do not mix the creation 
and annihilation operators a\ and a*, therefore they preserve Jq and are the generators of 
Isot{Jo). 

Transformations generated by Ea and e}- 

J ZJ 


f/ = exp( 


(2.46) 


necessarily mix creation and annihilation operators, U a U~^ = and provide a 

representation of the Squeeze{Jo) symplectic transformations. The corresponding Bogolyubov 
coefficients are 

sinh 


= cosh \/ /3/3t 


T = 




/ 3 . 


(2.47) 


A convenient parametrization of Squeeze (Jo) is provided by the Siegel domain S of com¬ 
plex symmetric matrices 7 in the unit disk [35], i.e. 

5 = {7 G Mat(A', C)| 7 = 7 * and 1 — 77 '^ > 0 } . (2.48) 

The squeeze unitary transformation (2.46) can be put in the form 


C/=exp( 


V 


- 7 *^' F^. 

2 ^ 


= exp 

where we have defined the matrix 


exp 


log(l - 77 ^)*^ Eij^ exp ^ 7 ”^ Fij^ 


tanh-y/^ 


(2.49) 

(2.50) 

(2.51) 


Note that 7 is symmetric and satisfies 0 < 77 ^ < 1 because of the properties of the hyperbolic 
tangent, therefore it belongs to the Siegel disk S. We can write the Bogoliubov coefficients 
in terms of 7 , 

= (1 — 77 l )“2 j T = (1 — 77'^)“2 q'. (2.52) 

Note that 7 = Using the expression for the symplectic transformation in terms of the 

Bogolyubov coefficients. 


My = '^{^ijUiU^j + '^ijUiuf + + ^*ijU*uf) 


(2.53) 




we can also write the symplectic matrix M in terms of 7 . By definition the reference complex 
structure Jq transforms non-trivially under conjugation by it is squeezed. We label a 
complex structure with a matrix 7 and define 

r—1 


Jo m: 


(2.54) 
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The space of complex structures can therefore be described by a reference complex 

structure Jq and a point 7 in the Siegel disk S. 


3 Squeezed vacuum states: properties 

Squeezed vacuum states can be characterized in a variety of equivalent ways. Here we follow 
the approach of [9] and show that squeezed vacua are group-theoretical coherent states for 
the symplectic group. In particular we show that they provide a resolution of the identity 
for the even sector of the Hilbert space of the system. Moreover we discuss the squeezing of 
their correlations. 


3.1 Vacuum state associated to a complex structure 

A squeezed vacuum associated to the complex structure J ^ J (M, fl) is a state in the Hilbert 
space Hv that satisfies the defining condition 

^^eiJ) = o. (3.1) 

Analogously to what we have done for the reference complex structure Jq we can introduce 
a set Ui of orthonormal eigenvectors with eigenvalue -|-i, i.e. Jui = +iui. The associated 
annihilation operator is a* = and the squeezed vacuum satisfies aj|J) = 0. Using the 

parametrization in terms of a point 7 in the Siegel disk we can also write 

(Oi-b 7 ijoJ) |J..y) = 0, (3.2) 

where a, and a\ are annihilation and creation operators associated to the reference complex 
structure Jq. Using Eq. (2.35) we see that the state \J^) can be obtained by squeezing the 
vacuum | Jq) with a unitary transformation U{M^), 


\J^) = UiM^)\Jo). (3.3) 

Using Eq. (2.50) together with EijIJq) = ^Sij\Jo) and Fij | Jq) = 0, we find the expression 

iXy) = det{l -exp \Jo) . (3.4) 

This expression allows us to compute the scalar product of two squeezed vacua 

det(l - 717 !) det(l - 7272 )^ ^ 


( J.yi I Tyj) ~ 


det(l - 7 | 72)2 


(3.5) 


In particular (Jo|Ty) = det(l — 77 ^) 4 . Moreover we can expand a squeezed vacuum over the 
basis (2.30), 

I Ty) = ( 7 ) |nj; Jq) , (3-6) 
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with coefficients 


Cmil) 


det(l — 77 '!') 4 


TT 


N 


[ 

Jc^ 




N 

n 

i=l 



(3.7) 


Note that the coefficients with odd all vanish: squeezed vacuum states 

have support on the sector of the Hilbert space spanned by Fock space basis elements In*; Jq) 
with even. 


3.2 Squeezed vacua as symplectic coherent states 

Group-theoretical coherent states are defined in terms of a Lie group G and its unitary 
irreducible representation U{g) on a Hilbert space [9]. Given a reference state |no) in 
a set of states is obtained by applying a group transformation, \vg) = U{g)\vQ). A set of 
coherent states is obtained by an appropriate choice of the reference state. 

Let g be the algebra of the group G, gc its complex hull and Xk the set of generators of 
gc. A generator preserves or annihilates the state |?;o) if there exists a Afc G C such that 

= Afc |no) . (3.8) 

Within the Hilbert space % there is a maximal set {Xi ,..., A„} of generators such that at 
least one solution of Eq. (3.8) exists. The reference state |no) for the construction of coherent 
states is required to be one of these solutions.® A case of particular interest occurs when 
a subset of the generators forms a subalgebra f} of gc. In this case the reference state 
satisfies 

U{g)\vo) = e^^\vo) \f g € Isot{ho). (3.9) 

i.e. it is invariant up to a phase under the action of the isotropy group Isot{ho) = {g G 
G\ghog~^ = ho} for some group element ho- In this case a coherent state is labeled by a 
reference state \vo) and an element of the coset G/Isot{ho)- This construction can be imme¬ 
diately applied to the unitary representation of Sp{2N, M) discussed in the previous section 
to obtain symplectic coherent states of a bosonic lattice. 


The Hilbert space Tiv of our bosonic system provides a reducible unitary representation 
of the double-cover of the symplectic group Sp{2N,M.). Defining the parity operator H = 
exppvri ~ 5 )] decompose T-Ly in 

T-Lv — T^even ©Hodd- (3.10) 

where Tieven is the space of eigenstates of H with eigenvalue -|-1 and 'Hodd with eigenvalue —1.^ 
As the unitary representation U{M) = exp (—z is constructed out of quadratic op¬ 

erators it leaves the parity of the state invariant and T-Ly is therefore reducible. An irreducible 
representation is induced on the even and the odd sector of the Hilbert space. 

®This condition guaranties that the set of coherent states have minimal dispersions of the Casimir operators 
of the group [9] . 

^The sector Heven is spanned by the basis elements (2.30) with ni + • • • -I- un even, Hodd by ni -I- • • • -I- njv 
odd. 
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Within Tieven there is a unique state | Jq) that satisfies the condition (3.8) for the maximal 
set of generators Eij and Fij, 

^ij\Jo) = ^^ij\Jo) , and Fij | Jq) = 0. (3.11) 

This reference state is invariant up to a phase under the action of the isotropy group of the 
element element Jq of Sp{2N,R). A symplectic coherent state in 'Heven is obtained by acting 
on I Jo) with a symplectic transformation U{M), i.e. 

\Jm)^U{M)\Jo). (3.12) 

Because of the invariance of the state | Jq) under the action of Isot{Jo), we can label a sym¬ 
plectic coherent state by the state | Jq) and an element of the squeeze coset Squeeze{Jo) = 
Sp{2N,'R)/Isot{jQ). This is exactly the squeezed vacuum state labeled by the complex struc¬ 
ture J = MJqM~^. 

Within J^odd there is no state that satisfies Fij\vo) = Xij\vo) for all i,j. There is however 
a set of N orthogonal states that satisfy Fij\vo) = 0. Defining the single-excitation state 
|lfc; Jo) = al\Jo) we have 

Fij {Ik] go) = 0 , (3.13) 

with k = 1... A^. By choosing one of them as reference state, we have the set of symplectic 
coherent states in J^odd given by 

lU; Jm) = U{M) 4 I Jo) = 4 I Jm) , (3.14) 

where 4 = U{M)a\U{M)~^. 

A remarkable property of group-theoretical coherent states is that they provide a resolu¬ 
tion of the identity as in integral over the group. For symplectic coherent states we have 

1 = / |-1m)(«^m| dp{M) + f |lfc; Jm)(1a;; «^m| dp{M), (3.15) 

JSp{2N,R) JSp{2N,R) 

where dp,{M) is the Haar measure on Sp{2N,M). In Eq. (3.15) the first term projects on 
Heven) the second on Hodd- Symplectic coherent states provide an overcomplete basis of H 
and every state can be written as a linear superposition of | Jm) and |li; Jm) for instance. 
Clearly the integration over M in Eq. (3.15) is redundant: elements of the compact subgroup 
Isot{Jo) ~ U{N) leave the integrand \Jm){Jm\ invariant. We can parametrize the integral 


over the Squeeze{Jo) coset by using the Bloch-Messiah decomposition; we write M = RTq R' 
with R, R' in Isot{Jo) and To diagonal and in Squeeze{Jo). In particular we have 

H(To)=exp(^s,(4-E,,)) (3.16) 

i 

with squeezing parameters Sj > 0 and Tq = with Xq = Diag(si,..., sn, —si ,..., —sjy). 
Eollowing [39] we find the induced invariant measure 

dp{To) = (n4j=i I coshsj - coshsjj) Hidsi. (3.17) 
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The following expressions provide equivalent representations of the resolution of the identity 
in T^even- 


leven = [ \J){J\ dfl{J) = [ \Jm){Jm\ (3.18) 

Jj(n,V) JSp(2iV,R) 

= [ \JrTo){JrTo\ dfj,{R)dfi{To) (3.19) 

Ju{N)xR^ 

where d^{J) can be understood as the invariant measure on the space of positive compatible 
complex structures (or over the Siegel disk S) and d^{R) is the invariant measure on the 
unitary group U{N). 

3.3 Correlations and nncertainty relations 

The expectation value of the elementary operator on a squeezed vacuum |J) is always 
vanishing, 

(J|e“|J)=0. (3.20) 

This is an immediate consequence of the fact that ^“|J) belongs to T-LoAA and is therefore 
orthogonal to \ J). On the other hand the two-point correlation function is non-trivial and 
completely determined by the complex structure J, 

= {J\ee\J) = i . (3.21) 

This relation can be derived for instance by noticing that + ifl“^) where 

5“^ = and (JieeV) = M“,M^(Jo|eeVo)- 

Note that the set of vectors Va in the vector space V'' = satisfying the condition 
C^^Va Vb = 1 form an ellipsoid. In fact only the symmetric part of the correlation function 
determines the ellipsoid. This part coincides with half of the inverse metric g°‘^ = J°^c 
introduced in Eq. (2.8), so that the correlation ellipsoid is defined by 

= 1 (3.22) 

Using the Bloch-Messiah parametrization of J = MJqM~^ with M = RTq we find that the 
principal axis of the ellipsoid come in pairs given by the vectors R°^i and R°‘i+N that span a 
symplectic plane. The axis have length where Si are the squeeze parameters introduced 
in Eq. (3.16). The ellipsoid associated to the state | Jq) is a sphere of radius 1, i.e. Si = 0: 
squeezing the state | Jq) with a unitary transformation U{M) produces a squeezed vacuum 
I J) with ellipsoid obtained by squeezing the sphere with the same symplectic transformation 
M. 
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Consider a pair of observables q = Qa^'^ and p = pai°‘ in T-Ly that are canonically conju¬ 
gate,® i.e. [g,p] = i. The dispersion® of q and p on the squeezed vacuum | J) is 

Ag = y/i J“c qaQb , = \JI PaPb ■ (3.23) 

Using the fact that is a positive matrix of determinant 1 when restricted to a sym- 

plectic subspace we can show that 

Ag Ap = ^ y/f qaPbY ■ (3.24) 

In particular, the Heisenberg uncertainty relations for the canonical couple q,p are saturated 
only when J°'c qaPb = 0 and are squeezed otherwise. 

4 Entanglement spectrum and entropy of squeezed vacua 

We are interested in describing properties of a subsystem. In terms of the symplectic vector 
space {V, H), we consider a subspace A = of dimension 2Na- Its symplectic complement 
is defined in Eq. (2.3). We require that H is a symplectic vector space with respect to 
the induced symplectic structure, (H, H). This requirement corresponds to the condition 
A n = {0}. The vector space B = is also symplectic. Therefore we have a symplectic 
decomposition of V, 

V = A®B, (4.1) 

with B of dimension 2Nb and Na + Nb = N. Without loss of generality, we assume that A 
is the smallest of the two, Na < Nb- We denote Pa and Pb : U —>■ U the projectors from V 
to the subspaces A and B respectively. 

At the level of the Hilbert space Tiy the symplectic decomposition (4.1) corresponds to 
a tensor product structure 

nv = nA®nB. ( 4 . 2 ) 

In this section we derive a Schmidt decomposition of symplectic matrices analogous to 
the one discussed in [20] and write the reduced density matrix 

PA^Ttb{\J){J\) (4.3) 

of a squeezed vacuum in terms of the associated entanglement Hamiltonian Ha- The entan¬ 
glement entropy of the state | J) restricted to the subsystem A is the von Neumann entropy 
of the reduced density matrix, 


Sa{\J)) ^-TAA{pA'^ogpA), (4.4) 

®These observables correspond to a symplectic basis in a two-dimensional symplectic subspace of (UH), 
i.e. Q.°‘^qaPb = 1. 

®Tlie dispersion of an observable O on a state |s) is defined as AO = \/ (s|0^|s) — (s|0|s)(s|0|s). 
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where the trace is over the Hilbert space T-La- The entanglement entropy can be easily 
computed from the spectrum of the entanglement Hamiltonian [40]. Following [18] we show 
that, at given two-point correlation function, squeezed vacua are the states with the largest 
entanglement entropy, a special case of the more general result of [41]. The expression of 
the entanglement entropy in terms of the projection of the complex structure is derived in 
Sec. 4.4. 

4.1 Normal modes and decomposition of the complex structure 

Given a complex structure J on F we are interested in its projection to the symplectic sub¬ 
space A, i.e. Pa JPa- In general the projected complex structure does not induce the complex 
structure on the subspace. Here we use Williamson theorem [20, 34-36] to put the projection 
of a complex structure in canonical form. 

Williamson theorem states that a positive-definite symmetric matrix K G Sym_|_(2iV, M) 
can be put in a canonical diagonal form via a symplectic transformation, i.e. K = M (A©A)M* 
where A = Diag(Ai,..., Aat) with A* > 0 and M G Sp{2N,M.). The invariants A* are the 
normal modes of K. We provide a constructive proof of this theorem. 

Given K G Sym_|_(2N, M) consider the hermitian matrix H = '/K \Jq y/K. As = H 
and H* = —H we have that it has real eigenvalues with opposite sign associated to orthogonal 
eigenvectors that are complex conjugate, i.e. 

Hvi = XiVi and Hv* = —XiV* (4-5) 

with Aj > 0 and i = 1... N. The matrix H is diagonalized by a 2N x 2N unitary transfor¬ 
mation U = {vi ■ ■ ■ Vn vl ■ ■ ■ u*) so that 

FT = [/ (A © -A) U-^ (4.6) 

with A = Diag(Ai,..., Xm) and (A © —A) the block diagonal matrix with A and —A on the 
diagonal. Let us define the unitary matrix 



Note that, because of the structure of the eigenvectors of H, the matrix R = Wlf/I is real 
and orthogonal, RR^ = 1. We define now the matrix M 

M = VkUW . (4.8) 

It is immediate to check that MJqM^ = Jq, therefore M G Sp{2N,'R). This concludes the 
constructive proof of Williamson theorem: a positive-definite symmetric matrix K can be 
written in terms of the normal modes A* as 

K = M{k® K)M^ (4.9) 
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where the invariants A* are the absolute value of the eigenvalues of the hermitian matrix 
H = ^/K iJq^/K and the symplectic matrix M is given by Eq. (4.8). We note also that 
the matrix y/KH\/K~^ is diagonalizable and it has the same eigenvalues of the matrix H. 
Therefore the Williamson invariants A* can be directly computed as the absolute value of the 
eigenvalues of the matrix \KJq. 


A complex structure J G J{V,Vl) can be obtained from the reference complex structure 
Jo by conjugation with a symplectic transformation, J = MJqM~^. Using the polar decom¬ 
position (2.16) we can write M = TR and J = T'^Jq, where G Sym^(2A^, M). We adopt the 
notation [PaMPa] ; A —)• A to denote the full-rank matrix obtained from PaMPa ■ V —)> 

By construction, the projection PaJqPa of the reference complex structure Jo to the subspace 
A is also a complex structure that we will still call Jo = [PaJqJA]. The projection of J on 
the other hand does not belong in general to J (A, Q) and can be written as 

[PaJPa] = [PaT‘^Pa]Jo. (4.10) 

with \Pa T‘^Pa\ G Sym^(2A^yi, M). We apply Williamson theorem to this matrix and we obtain 

[Pa T^Pa] = Ma (A © A) M\ (4.11) 

where Ma G Sp(2iVA,M) and A = Diag(cosh2ri,..., cosh2rjv^). From Williamson theorem 
we know that the normal modes cosh 2 rj can be determined as the absolute value of the 
eigenvalues of [PaT‘^Pa\ Uq- This is exactly i times the projection of the complex structure, 
so that we have [Pa UTA] = Ma (A © —A) i.e. 

Eigenvalues ([Pa i JPa]) = {=t cosh 2 rj} . (4.12) 

It follows that 

Eigenvalues^[P a ^ ^ Pa]J ={(coshrj)^, —(sinhrj)^}. (4.13) 

The projection of a complex structure is a complex structure only if the “squeeze” parameters 
ri are all vanishing. 

The same analysis can be performed for the projection Pb iJPe to the subspace B. 
The projected complex structure has the same non-vanishing eigenvalues {±cosh2rj}, plus 
2{Nb — Na) eigenvalues ±1. In particular, as for the projection to A, we can write 

[Pb UPb] = Ms (A © - A © 1 © - 1 ) M-^ (4.14) 

basis of E = A(B B consisting of the union of a basis of A and a basis of B can be chosen. In this case 
the matrix elements of [PaMPa] coincide with the matrix elements of M on a restricted range. We adopt this 
more abstract notation because it simplifies some algebraic manipulations when we use an orthonormal basis 
that is not adapted to the decomposition A® B. 
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where Mb G Sp(2A/^B,M). When extended to V we find that the symplectic transformation 
Ma © Mb G Sp(2A^, M) puts the complex structure in the canonical form 


J = {Ma® Mb) Tq Jq {M^^ © 
with To a symplectic transformation in Squeeze{Jo) given by 


Na 


i=l 


with 


1 

0 

0 ' 

II 

Ki 

0! 0 

1 


0 

0 


To = exp ( ^ r-j KiJo 
( 0 


Kj = 


V 


or more explicitly 



/ c © c 

s © —s 0 

To = 

s © —s 

c © c 0 



0 1 ^ 


0 / 


with c = Diag(coshri,... ,coshrAr^) and s = Diag(sinhri,... ,sinhrAr^). 


4.2 Schmidt decomposition and the entanglement Hamiltonian 

Using the decomposition 

M = {Ma © Mb) Tq R 

we can write the squeezed vacuum state associated to the complex structure J = 
in the form 


1 ^) 


U{M)\Jo) = U{Ma® Mb)U{To)\Jo) 


Na 

UaUb ex.p {'^ri{albl - aibi))\Jo) 

i=l 


Na 

uaub n 

i=l 



(tanhn)^ (4)^ {b\)^ 
cosh ri 


\Jo) 


where bi = aNA+i- Defining 


Na Na+Nb 

\JA) = U{MA)(^\jo)i , \Jb) = U{Mb) (g) Uo)* 

i=l i=NA+i 


and Ui = Ua di Ua^^ bi = UBbi we have the orthonormal basis in Ra and Rb 


^ (air 




|ni,.. .^una'^Ja) = ( H ^ ) I’^a) , |ni, • • • -^5)2 = ( H 




i=l 


\Jb) 


(4.15) 

(4.16) 

(4.17) 

(4.18) 

(4.19) 
MJoM-i 

(4.20) 

(4.21) 

(4.22) 

(4.23) 

(4.24) 
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This is a Schmidt basis for the squeezed vacuum | J), i.e. 


1 ^) 


E 

ni -- n ,^^=0 



(tauhri)*^* 
cosh r j 


ni,..., Ja) ® |ni,..., Jb) ■ 


(4.25) 


The reduced density matrix pA can now be easily computed, 


PA = Tib 


oo / Na 

E n 


ni'--njv^=0 \i=l 

This expression can also be written as 


(tanh ri) 
cosh r,; 


\ni,..., riAT^; JA){ni,..., Ja\ ■ 


PA = 


where Ha is the entanglement Hamiltonian 


Na 


^-2-kHa 

Za 


(4.26) 

(4.27) 


Ha =-log(tanh 

vr 


Ti) ajoi 


(4.28) 


i=l 


and Za = Tr(e“^-'^) = ])([^(coshrj)^. The squeezing parameters r* are obtained from the 
eigenvalues ±cosh2rj of iJ restricted to region H, Eq. (4.12). The entanglement spectrum is 
the set of eigenvectors and eigenvalues of the entanglement Hamiltonian as computed above 
[40]. 


4.3 Entanglement entropy and extremality of squeezed vacua 

The entanglement entropy of the squeezed vacuum | J) restricted to the sub Hilbert space Ha 
is given by the von Neumann entropy of the reduced density matrix, 

5A(|J)) = -Tr(pAlog/9A). (4.29) 

Let us consider two squeezed vacua | Ji) and | J 2 ) on Hy. Eor a given subsystem A we want 
to compare the entanglement entropy of the two states. It is useful to define the correlation 
matrices Cf^ and of the two states using Eq. (3.21), together with the entanglement 
Hamiltonians appearing in the expression of the reduced density matrices pi, p 2 of the two 
states. 

Hi = iKiabCC’^ H 2 = \K2abei" • (4.30) 

The relative entanglement entropy Ra is the relative entropy of the reduced density matrices 

.Ry4(|^2); jjl)) = Tr(^/92log/92 - /02logpi^ >0. (4.31) 
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It is always positive if the two density matrices do not coincide and zero if they do [18]. The 
difference between the entanglement entropies can be bounded from below using the following 
inequality 

5a(|Ji))-5a(|J 2)) =-Tv(pilogpi) +Tr(p 2 logpi) +i?A(|J 2 );|Ji)) (4.32) 

>-Tr(^pilogjOi^+Tr(^P2logpi^ (4.33) 

= Tr{H, pi) - p 2 ) (4.34) 

= p,) - TiiC^%2)) (4.35) 

= lK,ab{Cf-Cf) (4.36) 

Therefore a larger correlation function, > 0, implies a larger entropy for every 

subsystem A. In terms of the complex structures we have 

(Ji->0 ^ 5 a(|Ji)) >5a(|J2)). (4.37) 

The same method can be used to prove that squeezed vacuum states are the states with 
the largest entanglement entropy at fixed correlation function [41]. Consider a generic state 
|s) G Tiv and define the connected correlation function of this state in the standard way 

cf = {s\ee\s) - {s\as){s\e\s). ( 4 . 38 ) 

Out of we can define a complex structure 

= (4.39) 

The associated squeezed vacuum | J) that has exactly the same correlation function C'f' as 
the state |s). Then we have 

S’a(|J)) - S'a(|s)) = - Tr(pjlogpj) + Tr(pslogpj) + i?A(|s); |J)) (4.40) 

= ^KablrriCe PJ) Ps)) + RaUsYAJ)) (4.41) 

= \K^b{Cf - Cf) + RaHs)-, \J)) (4.42) 

= RAi\s);\J))>0 (4.43) 

Therefore at fixed correlation function the squeezed vacuum is the state with the largest 
entanglement entropy. 

4.4 Entanglement entropy of squeezed vacua in terms of the complex structure 

The entanglement entropy of the squeezed vacuum \J) restricted to the sub Hilbert space 
Ha is given by S'a(|J)) = — Tr(/3Alog/0A)- Using Eq. (4.27) for the reduced density matrix 
expressed in terms of the entanglement Hamiltonian Ha it is immediate to show that 

Na 

- TrA {pa log pa) = E ((coshrj)^ log(coshrj)^ — (sinhrj)^ log(sinhrj)^) . (4.44) 
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This expression can be cast in terms of the eigenvalues of Pa^-^Pa, Eq. (4.13). Equivalently, 
in matrix form, we find our final expression for the entanglement entropy in terms of the 
complex structure: 


Sa{\J)) = It^Pa^^i^Pa log 


Pa- 


1 — i J 


-Pa 


where the trace is over 2N dimensional matrices and not on the Hilbert space Pa- 
the main formula studied in the rest of the paper. 


(4.45) 
This is 


Note that this formula is equivalent to the one found in [19] for the entropy of a Gaussian 
bosonic channel expressed in terms of its symplectic invariants, and the one of [42] for the 
entanglement entropy of a quasi-free state of a scalar field on a causal set, expressed in terms 
of the covariant correlation function. The expression in terms of the complex structure makes 
manifest the fact the entanglement entropy vanishes when Pa^P^Pa is a projector (in this 
case in fact its eigenvalues are either 0 or 1 and OlogO = 0, llogl = 0). This happens if 
and only if the projection of the complex structure Pa JPa induces a complex structure for 
the subsystem A. To have a nonvanishing entanglement entropy we need that the complex 
structure J sends a vector in A into a vector with components in H. In fact, using Pa-\-Pb = 1 
and {Pa'^JPa)'^ = Pa + PaJPbJPa, we see that the eigenvalues of Pa^JPa are different 
from ±1,0 if and only if Pb JPa is non-vanishing. 


4.5 Entanglement entropy for small and large squeezing 

Suppose that the squeezing parameters r* in Eq. (4.12) are all small. This happens if Pa JPa 
is close to being a complex structure, in particular tr((P^iJP^)^ — l) <C Na- Expanding 
Eq. (4.44) r* <C 1 we find 


Na 

Sa{\J)) = - ^ogrf)r‘f ± 0{rf). (4.46) 

i=l 

On the other hand for large squeezing r* ^ 1, i.e. ti(^{P a^JPa)‘^ — l) 3> Na, we find the 
asymptotic expansion 

N 

Sa{\J))- [Y.^rij+{l-2log2)NA. (4.47) 

i=l 

The leading order can be written in matrix form in terms of the projected complex structure 
as 

<5a(|J)) -^log\deti[PAiJPA ])\, (4.48) 

where [Pa ^JPa] is the full rank matrix obtained from the projected symplectic structure. 
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5 Scalar field on a lattice and squeezed vacua 


In this section we apply the formalism of squeezed vacua and complex structures to study the 
entanglement entropy of a scalar field on a lattice. The lattice provides a notion of locality 
for the interaction between the N bosonic degrees of freedom. 

The entanglement entropy of the ground state of a local Hamiltonian [10-12] is well- 
studied [14-17] and know to lead to an area law [13]. In Sec. 5.2 we consider a generic lattice 
and use formula Eq. (4.45) to determine the dependence of the area law on the speed of sound, 
in the limit of small speed of sound. 

The time dependence of the entanglement entropy for a state that is not an eigenstate of 
the Hamiltonian has been studied using a variety of methods [43-48] . For a chaotic system the 
rate of the entanglement entropy growth is expected to be related to the Lyapunov exponents 
of the system [49-53]. In Sec. 5.3 we study the evolution of a squeezed vacuum in the presence 
of an instability described by a quadratic Hamiltonian with potential that is unbounded from 
below. We find a precise relation between the long-time behavior of the entanglement entropy, 
the size of the subsystem, the eigenvectors of the interaction potential and the eigenvalues 
associated to the unstable directions. In particular we find that the rate of growth of the 
entanglement entropy is asymptotically bounded from above by the Kolmogorov-Sinai rate 
of the associated classical system. 


5.1 Bosonic lattice: degrees of freedom and regions 


Consider a system of N bosonic degrees of freedom (</>j,7ri) with canonical commutation 
relations 

[(l)i,7rj] =i6ij , [(t>i,(l)j] = 0 , [Tri,TTj]=0, i,j = l,...,N. (5.1) 

The canonical variables {(j)i,TTi) are defined by 4>i = and tt* = 5i+]\fa^°' in terms of the 
observable on 77y. We can define a set of annihilation and creation operators as 


(*/’i T i VTj) , 


4 = 


(5.2) 


These are the operators associated to the reference complex structure Jq = ©^iJo introduced 
in Eq. (2.10). We denote | Jq) the vacuum annihilated by Oj and [uj; Jq) the Fock basis built 
over this vacuum as in Eq. (2.30). 

In order to define a notion of locality and adjacency relations between the degrees of 
freedom we equip the system with a lattice structure. We introduce a graph T with N nodes: 
each node plays the role of site i where the degree of freedom {(j)i,TTi) lives. 


Consider the set Af = (1,2,... ,N). We call nodes the elements i £ Af and we equip 
functions Vi on this set with the structure of a vector space. The adjacency relation of two 
nodes i,j is determined by the adjacency matrix Uj [54], a linear map from N io N defined 
by 


iij — 


1 

0 


ifl,j distinct and adjacent, 
otherwise. 


(5.3) 


We call link the unordered couple = (i,j) with i,j adjacent nodes. We denote C = 
{{hj)\h3 e = 1} the set of links. The couple T = {M,C) defines a finite graph that 

is undirected and simple. The number of links of the graph is determined by the adjacency 
matrix via the formula 

L = - ^ tij . (5.4) 

The valence di of the node i is defined as the number of links that contain the node i and is 
given by 

di = '^ tij . (5.5) 

isA/" 

The adjacency matrix provides us also with a notion of graph Laplacian defined as 

^ij — I'ij di 6ij . (5.6) 

The matrix — Ajj is positive semi-definite and its spectrum captures various properties of 
the graph. In particular in the following we assume that the graph T is connected^^. In 
this case there is at most one eigenvalue 0 in the spectrum of the graph Laplacian with 
the corresponding eigenvector being the constant vector. If we impose boundary conditions 
on the eigenvectors, such as vanishing components on a set of nodes of the graph, then all 
eigenvalues are strictly positive. 

We introduce now the notion of subgraph T^. Consider a subset A/a C Af with Na nodes. 
We denote Mb its complement so that M = A/a U Mb ■ We denote 

hidkj (5.7) 

the projector from M to A/a. We also define the subset Ca = £ A/a, tij = 1} C £ 

defined as the set of links connecting two nodes in A/a, and the set of boundary links 

k^dA = {(bi)K e j G Mb, tij = 1} . (5.8) 

We have £ = £a U £b U Cqa- The subgraph Ta C T is defined as Ta = (A/a, £a)- In general 
the subgraph Ta can have more than one connected component. 

A bosonic lattice (T, T-iv) is defined by a graph T with a bosonic degree of freedom {cpi, TTj) 
at each node i G M and the associated Hilbert space 

Hv = (g)£2(M). (5.9) 

ieJV 

A region of the bosonic lattice is a coucple {Ta,'Ha) where Ta is a connected subgraph of T 
with the bosonic degree of freedom ((/)j, vTj) with i G Ma- The Hilbert space of the region A is 

'Ha= (g) £'(M). (5.10) 

iGA/a 

graph F is connected if for every two nodes ii, *m there is a sequence [ii, 12 ,... Am) such that 
~ 
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Notice that we have the tensor product structure %v = T~La ®T-Lb- This decomposition 
corresponds to a decomposition of symplectic vector spaces V = A® B, where A = is 

the vector space spanned by the canonical basis {ej, JoCi} with Cj = (<5“) and z = 1 ,..., N^. 
The projector from y to A is 


Na 

Pa = ^ (eiej + (Joei)(Joeiyj . 


(5.11) 


2=1 


In terms of the projector Ha to the subset of sites Ma, Eq. (5.7), we have Pa = © II/i. 


5.2 Ground state of Hamiltonians and the area law 


Consider the lattice Hamiltonian 

H = ^ - (SijTTiTTj + Vij(l)i(t)j^ (5.12) 

with Vij(f>i(f)j a potential term bounded from below, i.e. V = {Vij) G Sym_,_(A/',M). The 
Hamiltonian can be put in normal form 

H = '^u)k{hlbk +^ (5.13) 

k=l 


where ujk are the eigenvalues of the matrix \/V and bk = Ukai°‘ = i is 

the annihilation operator associated to the complex structure 


J = 


0 

ve 0 y ■ 


(5.14) 


This expression can be obtained by writing the Hamiltonian in the form H = and 

applying Williamson theorem to Kab- The ground state of the Hamiltonian H is the squeezed 
vacuum | J) that satisfies bk\J) = 0. 

The entanglement entropy S'a(IT)) can be obtained by computing the eigenvalues A = 
±cosh2rj of the projected complex structure Pa ^JPa- The eigenvalues of Pa^JPa satisfy 


0 


det [Pa iJPa — a l) 


det 


/ -A1 

[iUAWUA 


— I 


UaW^Ua 

-XI 


= det (a^i- nA\/EnA\/i^nA) 

= det(nAVEns\/i^nA-(i-A^ji), (5.15) 

where we have used Ha + H^ = 1 in the last line. Therefore we find 

(sinh2rj)^ = —Eigenvalues^H aa/EHsa/E”^ Ha^ , 


(5.16) 







The squeezing parameters r* are non-vanishing and the entanglement entropy non-trivial only 
if the potential Vij entangles A with B, i.e. only if 7 ^ 0. 

In the following we consider a special case: the lattice discretization of the scalar field 
Hamiltonian on a trivial backgronnd, i.e. the potential term 

Vij = 6ij - c^Aij . (5.17) 

This potential corresponds to setting the lattice spacing and the mass of the scalar held to 
unity, while keeping the speed of sound c free. Note that the lattice is not assumed to be 
regular. It is instructive to study the behavior of the entanglement entropy for small speed 
of sound. For c <C 1 we have 

4 4 

- UaWUbVv^Ua = jVaA UbA Ua + 0(c®) = jT + 0(c®) (5.18) 

where T = HyitHstHA is the adjacency matrix tij restricted to the bonndary links Cqa- As 
a resnlt, using Eq. (5.16), 

4 

{rf} = Eigenvalues+ 0(c®). (5.19) 

Substituting in Eq. (4.46) we hnd 

Sn(|J))=tr((l-log|^T)|^T) +0{c^) (5.20) 

4-1 41 

= ^^log-^trT +0{c^) = ^ ^logArea(5H) + O(c^) (5-21) 

where tr T = ^^Cqa is the number of links of the graph that cross the boundary the region A 
i.e., in a lattice discretization, the area of the bonndary of the region A in nnits of the lattice 
cut-off. 

5.3 Linear generation of entanglement at instabilities 

We consider again a lattice Hamiltonian of the form 

H = '^ -(SijTTiTTj + Vij(j)i(j)j^ (5.22) 

but this time we do not assume that the Hamiltonian is bounded from below: the potential 
Vij(pi4>j term is assumed to have Nj unstable directions. We write 

Ni N 

V = ^-Xl Vkvi + ^ cjI Vkvi (5.23) 

k=l k=Nj-\-l 
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where are eigenvectors of the potential, — A| negative eigenvalues corresponding to the 
unstable directions and positive eigenvalues corresponding to the stable directions. We 

order the negative eigenvalues so that Ai > A 2 > • • • > Xnj > 0. We assume also that the 
eigenvalues are non-degenerate. We are interested in the time evolution of the entanglement 
entropy of an initial state | Jo). By writing the Hamiltonian in the form H = we 

find that the time evolution operator 

e-iHt ^ ^^ 24 ) 

corresponds to the symplectic transformation M = ex.p{—QKt) where QK = 

Therefore the evolution of the initial squeezed vacuum | Jq) is still a squeezed vacuum 

|Ji) = e-'^Vo) (5.25) 

with 

Jt = (5.26) 

This expression in principle allows us to compute the entanglement entropy 5^4 (|Jt)) as a 

function of time. Here we study the long-time limit under the assumption that there is at 
least one unstable direction, Ai > 0. 


The spectrum of Oit' can be expressed in terms of the eigenvalues and eigenvectors of 
the potential V. The eigenvalues /j of fiK satisfy the equation^^ 


0 = det(nK - fxl) = det(H + fiH), 


(5.27) 


therefore QK has 2N/ real eigenvalues /u = ±Afc and 2(N — Nj) imaginary eigenvalues n = 
Ticufc. The associated eigenvectors are Vk) and ■^{'^iijJkVk, Vk)- Using these results 

we can diagonalize the symplectic transformation M = exp(—QKt), i.e. 




M = NDN-^ 

(5.28) 

with D = Diag(e''"'^i*, e ^ 1 *,... 

p+A^Vt^ 

g+i'^Nj + lt g— 

, g-ioiivi) and^^ 

II 

toll 

' -AlUi 

-l-AiUi ••• -ilvnVn +iwArt'iV^ 

1 . (5.30) 


Vl ■■■ VN VN j 



19 /I 0\ f 0 V 

have K = { and fiA' = 

\0 V J \-l 0 

^^The first determinant is over 2A'^-dimensional matrices while the second on A^-dimensional matrices 
^■^The inverse of N is 


< -Ai ^vl v\ ^ 

-fA]^ ^v\ Vi 



+iuij/v%^ 

V*N 

■ —1 t 

t 


'^N / 


(5.29) 
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Note that 


N 

MM^ = Y,CkPk (5.31) 

k=l 

where = v^v^. is the projector on the eigenspace of the potential, is the 2x2 matrix 
(with detCfc = 1) 


Ck 


1 / 1 - A| + (1 + A|) cosh(2Afct) 

-(Afc + A^^)sinh(2Afct) 1 


-(Afc + A^^)sinh(2Afct) \ 

A^ ^ + (1 + A^ cosh(2Afct) / 


(5.32) 


and we have adopted the notation Afc = iw^ for /c = A^/ + 1 ,..., A^. 


The entanglement entropy S'a(| Jt)) depends on the projected complex structure Pa JiPa, 
Eq. (4.45). Under the assumption of large squeezing rj S> 1 (to be checked a posteriori), we 
can use the asymptotic expression (4.48), 

SA{\Jt)) ~ ^log|det([PAiJt^’A])| . (5.33) 

In the following we also assume, without loss of generality, that A is the smallest of the two 
subsystems in the bipartition, i.e. Na < Nb- Note that, using Jt = MJqM~^ = MM^Jq and 
[PaJqPa] = Jo, we have 

|det([PAiJtJAl])| = |det([PAMM*PA] iJo)| = |det([PAMM*PA])| . (5.34) 

Furthermore, using Eq. (5.31) and the spatial projector IIa (Eq. (5.7)), we obtain 

N N 

det([PA UiPa]) = det ^ HAPkP^A = det ^ , (5.35) 

k=l k=l 

where we defined Ck = CkP^APkP-A- We are interested in the asymptotic behavior of the 
entanglement entropy for large t, which requires the determination of the leading contribution 
to (5.35) in this limit. We write the entries of the matrices in the summand in (5.35) explicitly 
as 

[c,](-)JJ) = [CkT^ vivi , J = 1,..., ATa , a, /? = 1, 2. (5.36) 

Introducing an index r = 2(i — 1) + a, and noting that that the bi-index {i,a) is completely 
determined by r, we can switch to the single-index representation 

We recall now a standard result about the determinant of a sum of matrices as in 
Eq. (5.35). Using the definition of the determinant in terms of the Levi Civita tensors eri-.r 2 jv^ 
and the multilinearity of the determinant, we hnd 

/ AT \ N N 

det f j = ^ ••• ^ hki-k 2 N^)^ (5-37) 

\fc=l / k\=l fc2JV^=l 
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where bf^ki-k^N^) Siven by 


2Na 


2Na 


^{k-l---k2N a) — 


{2Na)\ 


E E 

n-r2iv^=lri-T-'j^^=l 


^n-raiv^er'-r'j^^ [CfeJ 




(5.38) 

In particular, for instance, 6(i...i) = detci. Let us focus now on a single term ^(A:r--fc2iVa) ™ 
Eq. (5.37). Note that, for fc = 1,..., N"/, each entry of the matrix Ck diverges as 

for large t, while the determinant det = 1 remains bounded. Since we are interested in 
the entropy for large t, it is sufficient to our purposes to keep track of the factors of 
in the calculation of the determinant. Now, the explicit form of the matrix given in 

(5.32) leads to several cancellations and simplifications in Eq. (5.38) for Consider 


for instance where there appear only two identical kiS. Naively one could expect 

an asymptotic behavior oc but it is in fact not the case. The reason is a cancellation 

that follows from det Cfc = 1. In the case of an odd number of identical k's appear, e.g. 
b{kkk■■■)■> the contribution is at most for large t. As a result, in the large t asymptotics 
of det([PA i-ATA]), each of the Nj exponential factors can appear at most once. 

In general, not all unstable modes k = 1,..., Ni contribute to the entropy. We define 


f 0 if n^PfcllA = 0 , r 0 if HBPk^B = 0 , 

bk ={ and hf = <^ (5.39) 

I 1 else, I 1 else. 

In particular 5^ = 0 corresponds to an eigenvector Vk that has its only non-vanishing com¬ 
ponents in B. In this case liAPk^A = 0, the corresponding k can be omitted from the 
sum in (5.35), and no factor of contributes to SA{\Jt)) in the asymptotic regime. On 
the other hand, if Vk has support in A only, i.e. 6^ = 0, we have a contribution that is 
not exponentially-diverging because of a cancellation of two terms that follows again from 
det Cfc = 1. Summarizing, a factor contributes to the asymptotics only if the eigenvector 
Vk associated to an unstable direction has non-vanishing components both in A and in B. 
Moreover each eigenvector with support in A only, i.e. 6^ = 0, reduces by 2 the total num¬ 
ber of factors that appear in the asymptotics. As a result, for large t, the asymptotic 
behavior is: 

N, 

|det([PAi^t-PA])| ~ exp Afc , (5.40) 

k=l 

where W is the smallest integer that solves 


Af* / Ni n \ 

= min 5^ 6^ , 2( hf) -2Nb\. (5.41) 

k=l V k=l k=l / 

The number of eigenvalues Afc contributing non-trivially to the asymptotics is 

and is always smaller or equal than the size 2Na of the subsystem. The eigenvalue Atv* 
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Figure 1. Entanglement entropy in the presence of instabilities. Left: entanglement entropy 5^(1 Jt)) 
as a function of time computed using Eq. (4.45). Right: time derivative of the entanglement entropy. 
The dashed line is the predicted asymptotic behavior from Eq. (5.43). The system considered consists 
of two degrees of freedom. The potential is + 4^2)'^■ We 

choose A ~ 1.5 for the unstable direction and uj = 1.1 for the stable direction. The subsystem is 
A = {(pi^TTi) and the initial state | Jq) is chosen to be the ground state of the potential Vij = Wq 5ij 
with frequency ujq = 1.0. The long-time behavior of the entanglement entropy is linear in time, with 
rate proportional to A and modulated with frequency w around the average. The asymptotic behavior 
is largely independent of the parameter wq of the initial state and provides numerical support for the 
relation to the Kolmogorov-Sinai rate hxs = A of Eq. (5.45). 


is the smallest element of the sequence Ai > A2 > • • • > Xni > 0 that contributes to the 
asymptotics. Using Eq. (5.33) we find that the long-time dependence of the entanglement 
entropy in the presence of instabilities is 

Af* 

SA{\Jt))-(Y.^t5^>^k) t. (5.42) 

k=l 


The subleading terms in the long time asymptotics depend on the oscillatory directions 
with normal frequencies ojk and result in the addition of a second term in Eq. (5.42). This 
term is bounded and oscillatory so that, defining the time-average S over the shortest period 
2ti jujk, we can write the asymptotics of the rate of entanglement entropy production as 


lim 

t^OO 


dSA{\Jt)) 

dt 


iv* 


k=l 


(5.43) 


Eig. 1 provides a concrete example of this behavior. Note that the classical system described 
by the Hamiltonian with potential Eq. (5.23) has local Lyapunov exponents given by ±Afc 
and Tiwfc [33]. The Kolmogorov-Sinai rate of the system is the sum over the real positive 
Lyapunov exponents [55], 

Ni 

hxs = > (5.44) 

fc=i 
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a quantity that is in general expected to measure the rate of growth of the coarse-grained 
entropy of the classical dynamical system [33, 55-57]. Formula (5.42) shows that the long¬ 
time behavior of the entanglement entropy in the presence of instabilities is bounded from 
above by the Kolmogorov-Sinai entropy rate, 


lim < hKS . (5.45) 

t —^OO (jt 

The equality between the entanglement-entropy rate and hxs is attained when each one of 
the unstable directions of the potential affects both the subsystem A and its complement 
B, and the number of unstable directions does not exceed the size of the bi-partition, Nj < 
min(2A^^, 2Nb)- 

6 Discussion 

We have shown that, treating squeezed vacua as coherent states for the symplectic group, a 
new formula for their entanglement entropy can be derived, Eq. (4.45). This formula involves 
only the complex structure J entering the definition of the squeezed vacuum and its projection 
Pa JPa to the phase space of the subsystem. We applied this formula to a scalar field on 
a generic lattice, investigated the relation between the area law and the speed of sound 
Eq. (5.21), and determined a precise relation between the entanglement entropy production 
in the presence of instabilities and the local Lyapunov exponents of the system, Eq. (5.43). 

The extension to a continuum quantum field theory of our results is non-trivial because 
of the presence of infinitely many degrees of freedom. The formulation in terms of complex 
structures is however well-suited for this analysis: the Fock spaces built on the squeezed 
vacua jJi) and IJ2) are unitarily equivalent if and only if tr(Ji — J 2 )“^ < 00 [8, 29]. In this 
case Eq. (4.45) is likely to result in a finite difference of entanglement entropies between the 
squeezed vacua jJi) and IJ2), therefore providing a generalization of the results derived in 
[30-32]. Moreover, formula (4.45) goes beyond the notion of geometric entropy [10-12]: the 
subsystem A does not have to be localized in space, it can be for instance a slice in momentum 
space [58]. In particular this remark applies to the analysis of the growth of entanglement 
entropy in the presence of instabilities of Sec. 5.3 and is directly relevant for the study of 
entanglement entropy production in inflationary cosmology [5-7]. 
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A Symplectic algebra and quadratic generators 


Sp{2N, M) is a connected Lie group. In this appendix we describe the algebra of its generators 
[34]. Writing M = expeA and expanding to the first order in e we find 

the relation = 0. This equation can be solved by posing X = KQ, with 

K = (Kab) a symmetric matrix. Therefore we have that M = exp(iLn) is a symplectic matrix 
and the symplectic algebra is 

sp(2A,M) = {AOj A e Sym(2A,M)}. (A.l) 

The commutator of two generators is 

[Kin, K2n] = KsQ (A.2) 

where = KinK 2 — K 2 nKi is also symmetric. The dimension of the symplectic algebra is 
dimsp(2A, M) = dim 5'?/m(2A, M) = N{2N + 1) . (A.3) 


The generators of the orthogonal group are skew-symmetric and therefore the matrices 
generating the Isot{Jo) subgroup have the form 

K,g = K- JoKJo (A.4) 


for some symmetric matrix K. This equation can be solved in block form 


Kis 


I Ai A2\ 
V-A2 Ai) 


(A.5) 


with Ai G Sym{N,M) and A 2 G Skew{N,M). The antisymmetric matrix KigJo generates the 
0(2iV, M) subgroup of the linear symplectic group. As a generic transformation in Isot{Jo) 
can be written as R = exp AisJo, we have 


dim/sot( Jo) = dim Sym{N,R) + dim Skew{N,R) = iV^ . (A.6) 


The symplectic matrices in Squeeze{Jo) are symmetric and positive, therefore the gener¬ 
ators have the form Ks^Jq with 

Kg^ = K + JoKJo (A.7) 

for some symmetric matrix K. This equation can be solved in block form 

with Bi, B 2 G Sym{N,'K). A generic transformation in Squeeze{Jo) can be written as T = 
expiLsqJo. Therefore we have 

dim Squeeze{jQ) = 2 dim Sym{N,R) = N{N + 1) . (A.9) 
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The exponential map from the Lie algebra to the group is not surjective. However, any 
element of the group can be generated by the group multiplication of two elements. Using 
the polar decomposition M = T Rwe can write any symplectic matrix M in a unique way as 
the product of a matrix in Squeeze{Jo) and a matrix in Isot{Jo), 

M = exp(iLsq exp(iLis 0). (A.10) 

This form will be especially important in the following sections where we discuss a unitary 
representation of the linear symplectic group. 

In terms of Poisson brackets, Sec. (2.3), symplectic matrices correspond to linear canon¬ 
ical transformations generated by a quadratic Hamiltonian H = In particular 

transformations in Isot{Jo) have generator 

PiPj +AYqiqj + {qiPj -piqj)^ (A.ll) 

with Ai G Sym{N,R) and A 2 G Skew{N,R) as in Eq. (A.5). Transformations in Squeeze{Jo) 
have generator 

CC’’ = ^ (^B\^ PiPj - qiQj + B^^ {qi pj + Pi Qj)^ (A. 12) 

with Bi,B 2 G Sym{N,R) as in Eq. (A.8). 

The analogous expressions in terms of creation and annihilation operators defined in 
Sec. (2.4-2.5) are 

Ei,, (A.IS) 

where a G Mat{N, C) is the hermitian matrix a = Ai -|-iA2 and /3 G Mat{N, C) the symmetric 
matrix /? = H2 — iHi. The Bogolyubov coefficients and Tjj can be expressed in terms of 
a and /3 as explained in Sec. (2.5). 
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